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1 Introduction 



Bogolubov's 1947 pairing theory for a Bose fluid was used by Foldy in 1961 to calculate the 



ground state energy of the one-component plasma (also known as "jellium") in the high density 
regime — which is the regime where the Bogolubov method was thought to be exact for this 
problem. Foldy's result will be verified rigorously in this paper; to our knowledge, this is the 
first example of such a verification of Bogolubov's theory in a three-dimensional system of bosonic 
particles. 

Bogolubov proposed his approximate theory of the Bose fluid [^] in an attempt to explain 
the properties of liquid Helium. His main contribution was the concept of pairing of particles with 
momenta k and —k; these pairs are supposed to be the basic constituents of the ground state (apart 
from the macroscopic fraction of particles in the "condensate" , or A; = state) and they are the 
basic unit of the elementary excitations of the system. The pairing concept was later generalized 
to fermions, in which case the pairing was between particles having opposite momenta and, at the 
same time, opposite spin. Unfortunately, this appealing concept about the boson ground state has 
neither been verified rigorously in a 3-dimensional example, nor has it been conclusively verified 
experimentally (but pairing has been verified experimentally for superconducing electrons). 

The simplest question that can be asked is the correctness of the prediction for the ground state 
energy (GSE). This, of course, can only be exact in a certain limit ~ the 'weak coupling' limit. 
In the case of the charged Bose gas, interacting via Coulomb forces, this corresponds to the high 
density limit. In gases with short range forces the weak coupling limit corresponds to low density 
instead. 

Our system has N bosonic particles with unit positive charge and coordinates Xj, and a uni- 
formly negatively charged 'background' in a large domain of volume V. We are interested in the 
thermodynamic limit. A physical realization of this model is supposed to be a uniform electron sea 
in a solid, which forms the background, while the moveable 'particles' are bosonic atomic nuclei. 
The particle number density is then p = N/V and this number is also the charge density of the 
background, thus ensuring charge neutrality. 

The Hamiltonian of the one-component plasma is 



where p = —iV is the momentum operator, = —A, and the three potential energies, particle- 
particle, particle-background and background-background, are given by 




(1) 



PP 




(2) 



l<i<jr<Ar 




(3) 
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Ubb = hp" / \x-y\-U''xd% . (4) 



2 

In Jn 



In our units /m = 1 and the charge is e = 1. The 'natural' energy unit we use is two Rydbergs, 
2Ry = me'^/JT?. It is customary to introduce the dimensionless quantity = {3/ Air p)^/^ e'^m/ It? . 
High density is small r^. 

The Coulomb potential is infinitely long-ranged and great care has to be taken because the 
finiteness of the energy per particle in the thermodynamic limit depends, ultimately, on delicate 
cancellations. The existence of the thermodynamic limit for a system of positive and negative 
particles, with the negative ones being fermions, was shown only in 1972 | LLe| (for the free energy, 
but the same proof works for the ground state energy). Oddly, the jellium case is technically a bit 



harder, and this was done in 1976 [LN| (for both bosons and fermions). One conclusion from this 
work is that neutrality (in the thermodynamic limit) will come about automatically - even if one 
does not assume it - provided one allows any excess charge to escape to infinity. In other words, 
given the background charge, the choice of a neutral number of particles has the lowest energy in 
the thermodynamic limit. A second point, as shown in [LN|, is that cq is independent of the shape 
of the domain $7 provided the boundary is not too wild. For Coulomb systems this is not trivial 
and for real magnetic systems it is not even generally true. We take advantage of this liberty and 
assume that our domain is a cube [0, L] x [0, L] x [0, L] with = V. 

We note the well-known fact that the lowest energy of H in (||) without any restriction about 
'statistics' (i.e., on the whole of (8)^L^(M^)) is the same as for bosons, i.e., on the symmetric 
subspace of (8)^L^(R^). The fact that bosons have the lowest energy comes from the Perron- 
Frobenius Theorem applied to — A. 

Foldy's calculation leads to the following theorem about the asymptotics of the energy for small 
rs, which we call Foldy's law. 

1.1. THEOREM (Foldy's Law). 

Let Eq denote the ground state energy, i.e., the bottom of the spectrum, of the Hamiltonian H 
acting in the Hilbert space Cgi^L^(R^). We assume that Q = [0, L] x [0, L] x [0, L]. The ground state 
energy per particle, cq = Eq/N , in the thermodynamic limit N,L ^ oo with N/V = p fixed, in 
units ofme^/h'^, is 

/A \ 1/4 

lim Eo/N = eo = -0.40154 r,-^/^ + o{p^^^) = -0.40154 ( — ) p^/^ + o{p^/^) , (5) 

>CXD \ 3 J 

where the number -0.40154 is, in fact, the integral 

A = r {p\p^ + 2)V2 - / - i| dp = - fXflill " -0-40154 . (6) 

TT Jq I ) 5V7rr(5/4) 

Actually, our proof gives a result that is more general than Theorem |1.1[ We allow the particle 
number N to be totally arbitrary, i.e., we do not require N = pV. Our lower bound is still given 
by (^), where now p refers to the background charge density. 
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In 0.40154 is replaced by 0.80307 since the energy unit there is 1 Ry. The main result of our 
paper is to prove (|5|) by obtaining a lower bound on Eq that agrees with the right side of (|5|) An 
upper bound to Eq that agrees with (|5|) (to leading order) was given in 1962 by Girardeau [ |GM |, 
using the variational method of himself and Arnowitt [ |GA| . Therefore, to verify (^) to leading 
order it is only necessary to construct a rigorous lower bound of this form and this will be done 
here. It has to be admitted, as explained below, that the problem that Foldy and Girardeau treat is 
slightly different from ours because of different boundary conditions and a concommitant different 
treatment of the background. We regard this difference as a technicality that should be cleared up 



one day, and do not hesitate to refer to the statement of 1.1 as a theorem. 

Before giving our proof, let us remark on a few historical and conceptual points. Some of the 
early history about the Bose gas, can be found in the lecture notes 0. 

Bogolubov's analysis starts by assuming periodic boundary condition on the big box and 
writing everything in momentum (i.e., Fourier) space. The values of the momentum, k are then 
discrete: k = (27r/L)(mi, m2, ^3) with nii an integer. A convenient tool for taking care of various 
n\ factors is to introduce second quantized operators (where a* denotes a or a*), but it has to 
be understood that this is only a bookkeeping device. Almost all authors worked in momentum 
space, but this is neither necessary nor necessarily the most convenient representation (given that 
the calculations are not rigorous). Indeed, Foldy 's result was reproduced by a calculation entirely 
in x— space |LS]. Periodic boundary conditions are not physical, but that was always chosen for 
convenience in momentum space. 

We shall instead let the particle move in the whole space, i.e., the operator H acts in the Hilbert 
space L'^(M.^^), or rather, since we consider bosons, in the the subspace consisting of the A^-fold 
fully symmetric tensor product of L^(M'^). The background potential defined in (|2|) is however still 
localized in the cube 17. We could also have confined the particles to Q with Dirichlet boundary 
conditions. This would only raise the ground state energy and thus, for the lower bound, our setup 
is more general. 

There is, however, a technical point that has to be considered when dealing with Coulomb 
forces. The background never appears in Foldy's calculation; he simply removes the k = mode 
from the Fourier transform, v of the Coulomb potential (which is i'{k) = 47r|A;|~^, but with k 
taking the discrete values mentioned above, so that we are thus dealing with a 'periodized' Coulomb 
potential). The k = elimination means that we set i^(0) = 0, and this amounts to a subtraction 
of the average value of the potential - which is supposed to be a substitute for the effect of a 
neutralizing background. It does not seem to be a trivial matter to prove that this is equivalent 
to having a background, but it surely can be done. Since we do not wish to overload this paper, 
we leave this demonstration to another day. In any case the answers agree (in the sense that our 
rigorous lower bound agrees with Foldy's answer), as we prove here. If one accepts the idea that 
setting z^(0) = is equivalent to having a neutralizing background, then the ground state energy 



problem is finished because Girardeau shows | GM | that Foldy's result is a true upper bound within 
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the context of the i/(0) = problem. 

The potential energy is quartic in the operators . In Bogolubov's analysis only terms in 
which there are four or two operators are retained. The operator Og creates, and oq destroys 
particles with momentum and such particles are the constituents of the 'condensate'. In general 
there are no terms with three Oq operators (by momentum conservation) and in Foldy's case there 
is also no four Og term (because of the subtraction just mentioned). 

For the usual short range potential there is a four term and this is supposed to give the 
leading term in the energy, namely cq = 47rpa, where a is the 'scattering length' of the two-body 
potential. Contrary to what would seem reasonable, this number, 4:iTpa is not the coefficient of the 
four term, and to to prove that Airpa is, indeed, correct took some time. It was done in 1998 



|LY] and the method employed in |LY] will play an essential role here. But it is important to be 
clear about the fact that the four Og , or 'mean field' term is absent in the jellium case by virtue 
of charge neutrality. The leading term in this case presumably comes from the two ag terms, and 
this is what we have to prove. For the short range case, on the other hand, it is already difficult 
enough to obtain the Anpa energy that going beyond this to the two Cg terms is beyond the reach 
of rigorous analysis at the moment. 

The Bogolubov ansatz presupposes the existence of Bose Einstein condensation (BEC). That 
is, most of the particles are in the k = mode and the few that are not come in pairs with 
momenta k and —k. Two things must be said about this. One is that the only case (known to 
us) in which one can verify the correctness of the Bogolubov picture at weak coupling is the one- 



dimensional delta- function gas [LLi] — in which case there is presumably no BEC (because of the 
low dimensionality). Nevertheless the Bogolubov picture remains correct at low density and the 
explanation of this seeming contradiction lies in the fact that BEC is not needed; what is really 
needed is a kind of condensation on a length scale that is long compared to relevant parameters, 



but which is fixed and need not be as large as the box length L. This was realized in |LY| and the 
main idea there was to decompose into fixed-size boxes of appropriate length and use Neumann 
boundary conditions on these boxes (which can only lower the energy, and which is fine since we 



want a lower bound). We shall make a similar decomposition here, but, unlike the case in |LY] 
where the potential is purely repulsive, we must deal here with the Coulomb potential and work 
hard to achieve the necessary cancellation. 

The only case in which BEC has been proved to exist is in the hard core lattice gas at half-filling 
(equivalent to the spin- 1/2 XY model) [ KLS| ]. 

Weak coupling is sometimes said to be a 'perturbation theory' regime, but this is not really so. 



In the one-dimensional case [LLi| the asymptotics near p = is extremely difficult to deduce from 
the exact solution because the 'perturbation' is singular. Nevertheless, the Bogolubov calculation 
gives it effortlessly, and this remains a mystery. 

One way to get an excessively negative lower bound to cq for jellium is to ignore the kinetic 
energy. One can then show easily (by an argument due to Onsager) that the potential energy alone 
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is bounded below by cq ~ —p^^^. See [LN|. Thus, our goal is to show that the kinetic energy 
raises the energy to — p^/^. This was done, in fact, in | CLY | , but without achieving the correct 
coefficient -0.803(47r/3)i/l Oddly, the -p^/^ law was proved in ]CLY[ | by first showing that the 
non-thermodynamic N'^^^ law for a two- component bosonic plasma, as conjectured by Dyson Q, 
is correct. 

The [CLY] paper contains an important innovation that will play a key role here. There, too, 
it was necessary to decompose into boxes, but a way had to be found to eliminate the Coulomb 
interaction between different boxes. This was accomplished by not fixing the location of the boxes 
but rather averaging over all possible locations of the boxes. This 'sliding localization' will play a 



key role here, too. This idea was expanded upon in | GG ]. Thus, we shall have to consider only one 
finite box with the particles and the background charge in it independent of the rest of the system. 
However, a price will have to be paid for this luxury, namely it will not be entirely obvious that 
the number of particles we want to place in each box is the same for all boxes, i.e., where i 
is the length of box. Local neutrality, in other words, cannot be taken for granted. The analogous 



problem in |LY] is easier because no attractive potentials are present there. We solve this problem 
by choosing the number, n, in each box to be the number that gives the lowest energy in the box. 
This turns out to be close to n = pi^, as we show and as we know from |LN] must be the case as 

Finally, let us remark on one bit of dimensional analysis that the reader should keep in mind. 
One should not conclude from that a typical particle has energy p^/^ and hence momentum 
Qj- Broglie wavelength p^^^^. This is not the correct picture. Rather, a glance at the 
Bogolubov/Foldy calculation shows that the momenta of importance are of order p~^^'^, and the 
seeming paradox is resolved by noting that the number of excited particles (i.e., those not in the 
A; = condensate) is of order Np~^^^. This means that we can, hopefully, localize particles to 
lengths as small as p~^^^^'', and cut off the Coulomb potential at similar lengths, without damage, 
provided we do not disturb the condensate particles. It is this clear separation of scales that enables 
our asymptotic analysis to succeed. 



2 Outline of the Proof 



The proof of our main theorem 1.1 is rather complicated and somewhat hard to penetrate, so we 



present the following outline to guide the reader. 



2.1 Section |3| 

Here we localize the system whose size is L into small boxes of size £ independent of L, but dependent 
on the intensive quantity p. Neumann boundary conditions for the Laplacian are used in order to 
ensure a lower bound to the energy. We always think of operators in terms of quadratic forms and 
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the Neumann Laplacian in a box Q is defined for all functions in ■0 G L'^iQ) by the quadratic form 



The lowest eigenfunction of the Neumann Laplacian is the constant function and this plays the role 
of the condensate state. This state not only minimizes the, kinetic energy, but it is also consistent 
with neutralizing the background and thereby minimizing the Coulomb energy. The particles not 
in the condensate will be called 'excited' particles. 



in the Introduction. The interaction among the boxes is controlled by using the sliding method 



background in each little box separately. 

The particles have to be distributed among the boxes in a way that minimizes the total energy. 
We can therefore not assume that each box is neutral. Instead of dealing with this distribution 
problem we do a simpler thing which is to choose the particle number in each little box so as to 
achieve the absolute minimum of the energy in that box. Since all boxes are equivalent this means 
that we take a common value n as the particle number in each box. The total particle number 
which is n times the number of boxes will not necessarily equal A^, but this is of no consequence 
for a lower bound. We shall show later, however, that it equality is nearly achieved, i.e., the the 
energy minimizing number n in each box is close to the value needed for neutrality. 

2.2 Section g 

It will be important for us to replace the Coulomb potential by a cutoff Coulomb potential. There 
will be a short distance cutoff of the singularity at a distance r and a large distance cutoff of the 
tail at a distance R, with r < R <^ i. One of the unusual features of our proof is that r are R are 
not fixed once and for all, but are readjusted each time new information is gained about the error 
bounds. 

In fact, already in Sect. ^ we give a simple preliminary bound on n by choosing R ~ p^^^^, 
which is much smaller than the relevant scale p~^^^, although the choice of R that we shall use at 
the end of the proof is of course much larger than p~^^^, but less than £. 

2.3 Section |5| 

There are several terms in the Hamiltonian. There is the kinetic energy, which is non-zero only 
for the excited particles. The potential energy, which is a quartic term in the language of second 
quantization, has various terms according to the number of times the constant function appears. 
Since we do not have periodic boundary conditions we will not have the usual simplification caused 
by conservation of momentum, and the potential energy will be correspondingly more complicated 
than the usual expression found in textbooks. 




Q 



To avoid localization errors we take i ^ p 



which is the relevant scale as we mentioned 



of [CLY|. The result is that we have to consider only interactions among the particles and the 
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In this section we give bounds on the different terms in the Hamiltonian and use these to get 
a first control on the condensation, i.e., a control on the number of particles n_|_ in each little box 
that are not in the condensate state. 

The difficult point is that is an operator that does not commute with the Hamiltonian and 
so it does not have a sharp value in the ground state. We give a simple preliminary bound on its 
average in the ground state by again choosing R ~ p^^l"^ . In order to control the condensation 
to an appropriate accuracy we shall eventually need not only a bound on the average, (n+), but 
also on the fluctuation, i.e, on (ra^). This will be done in Sect. ^ using a novel method developed 
in Appendix ^ for localizing off-diagonal matrices. 

2.4 Section § 

The part of the potential energy that is most important is the part that is quadratic in the con- 
densate operators Og and quadratic in the excited variables Op with p 7^ 0. This, together with 
the kinetic energy, which is also quadratic in the , is the part of the Hamiltonian that leads to 
Foldy's law. Although we have not yet managed to eliminate the non-quadratic part up to this 
point we study the main 'quadratic' part of the Hamiltonian. It is in this section that we essentially 
do Foldy's calculation. 

It is not trivial to diagonalize the quadratic form and thereby reproduce Foldy's answer because 
there is no momentum conservation. In particular there is no simple relation between the resolvent 
of the Neumann Laplacian and the Coulomb kernel. The former is defined relative to the box and 
the latter is defined relative to the whole of M'^. It is therefore necessary for us to localize the 
wavefunction in the little box away from the boundary. On such functions the boundary condition 
is of no importance and we can identify the kinetic energy with the Laplacian in all of "M?. This 
allows us to have a simple relation between the Coulomb term and the kinetic energy term since 
the Coulomb kernel is in fact the resolvent of the Laplacian in all of . 

When we cut off the wavefunction near the boundary we have to be very careful because we 
must not cut off the part corresponding to the particles in the condensate. To do so would give too 
large a localization energy. Rather, we cut off only functions with sufficiently large kinetic energy 
so that the localization energy is relatively small compared to the kinetic energy. The technical 
lemma needed for this is a double commutator inequality given in Appendix |^. 

2.5 Section 

At this point we have bounds available for the quadratic part (from Sect. ^) and the annoying 
non-quadratic part (from Sect. ^) of the Hamiltonian. These depend on r, i?, n, and (ra^). 

We avail ourselves of the bounds previously obtained for n and {n+) and now use our freedom to 
choose different values for r and R to bootstrap to the desired bounds on n and (n+), i.e., we prove 
that there is almost neutrality and almost condensation in each little box. 
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2.6 Section |8| 

In order to control (n^) we utilize, for the first time, the new method for localizing large matrices 
given in Appendix^. This method allows us to restrict to states with small fluctuations in n+, and 
thereby bound (n^), provided we know that the terms that do not commute with have suffciently 
small expectation values. We then give bounds on these n+ 'off-diagonal' terms. Unfortunately, 
these bounds are in terms of positive quantities coming from the Coulomb repulsion, but for which 
we actually do not have independent a-priori bounds. Normally, when proving a lower bound to a 
Hamiltonian, we can sometimes control error terms by absorbing them into positive terms in the 
Hamiltonian, which are then ignored. This may be done even when we do not have an a-priori 
bound on these positive terms. If we want to use Theorem A.l in Appendix we will need an 



absolute bound on the 'off-diagonal' terms and we can therefore not use the technique of absorbing 
them into the positive terms. The decision when to use the theorem in Appendix ^ or use the 
technique of absorption into positive terms is resolved in Section ^. 

2.7 Section § 

Since we do not have an a-priori bound on the positive Coulomb terms as described above we 
are faced with a dichotomy. If the positive terms are, indeed, so large that enough terms can be 
controlled by them we do not need to use the localization technique of Appendix |A| to finish the 
proof of Foldy's law. The second possibility is that the positive terms are bounded in which case 
we can use this fact to control the terms that do commute with n_|_ and this allows us to use the 
localization technique in Appendix |A| to finish the proof of Foldy's law. Thus, the actual magnitude 
of the positive repulsion terms is unimportant for the derivation of Foldy's law. 

3 Reduction to a small box 

As described in the previous sections we shall localize the problem into smaller cubes of size i <^ L. 
We shall in fact choose ^ as a function of p in such a way that p^^^i — > cxd as /) ^ oo. 

We shall localize the kinetic energy by using Neumann boundary conditions on the smaller 
boxes. 

We shall first, however, describe how we may control the electostatic interaction between the 



smaller boxes using the sliding technique of [CLY| 



Let t, with < t < 1/2, be a parameter which we shall choose later to depend on p in such a 
way that t — > as /? ^ oo. 

The choice of £ and t as functions of p will be made at the end of section ^ when we complete 
the proof of Foldy's law. 

Let X e C^(R3) satisfy supp X C [(-1 + 1)/2, (1 - t)/2f, < X < I, X{x) = 1 for x in the 
smaller box [(—1 + 2t)/2, (1 — 2t)/2]^, and x{x) = x{—x). Assume that all m-th order derivatives 
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of X are bounded by Cmt~"^, where the constants Cm depend only on m and are, in particular, 
independent of t. Let Xe{x) = x{x/t). Let r] = \Jl — X- We shall assume that X is defined such 
that r/ is also . Let = ri{x/l). Using x we define the constant 7 by 7^-*^ = J xivY dy, and 

note that 1 < 7 < (1 — 2t)^^. We also introduce the Yukawa potential Yi^{x) = \x\^'^e^'^^^^ for 
> 0. 

As a preliminary to the following lemma 3.1 we quote Lemma 2.1 in |CLY]. 
LEMMA Let K be given by 

K{z) = r-^ {e-'''' - e-^'h{z)] 

with r = \z\ and to > u > 0. Let h satisfy (i) h is a function of compact support; (ii) h{z) = 
1 A- ar'^ + O(r^) near z = 0. Let h{z) = h{—z), so that K has a real Fourier transform. Then there 
is a constant, C3 (depending on h) such that if u) — u > C3 then K has a positive Fourier transform 
and, moreover, 

eiCjK^Xi - Xj) > -{v - uj)N 

l<i<j<N 

for all xi, . . . xn G and all ei = ±1. 

3.1. LEMMA (Electrostatic decoupling of boxes using sliding). 

There exists a function of the form uj{t) = Ct^'^ (we assume that uj{t) > 1 for t < 1/2) and a 
constant 7 with 1 < 7 < (1 — 2t)^^ such that if we set 

w{x,y) = xdx)Y^(^tyf{x - y)Xe{y) (7) 
then the potential energy satisfies 



Upp + Upb + Ubb 

> 7 X] / ^'"{ Yl w{xi + {fi + X)i, Xj + (/i + \)€) 



A»tL o.oJ 



N 

P 

3 



V / w {xj + (/i + A)^, y + + X)£) dy 
.7=1 

j j w{x + {ijL + X)i, y + ifi + X)i) dx dy^ 



uj{t)N 

2e ' 



Proof. We calculate 



j di^ijxix+ip + \))YUx-y)x{y + {iJ. + X)) 
^^^Ve[-i/2,i/2p 



7X(x + z)Y^{x - y)x{y + z)dz = h{x - y)Y^{x - y), 
where we have set h = 7X * X- Note that h{0) = 1 and that h satisfies all the assumptions in 



Lemma 2.1 in [CLY|. We then conclude from Lemma 2.1 in |CLY] that the Fourier transform of 
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the function F{x) = \x\ ^ — /i(x)y^(f) (x) is non-negative, where a; is a function such that Lo{t) oo 
as t ^ 0. [The detailed bounds from [|CLY | show that we may in fact choose ui{t) = Ct~^, since 



uj{t) has to control the 4th derivative of h.] Note, moreover, that lim^^^o -^(2;) = uj{t). Hence 

N 



X — y) dxdy > 



The lemma follows by writing = F(x) + (x) and by rescaling from boxes of size 

1 to boxes of size i. □ 

As explained above we shall choose the parameters t and i as functions of p at the very end of 
the proof. We shall choose them in such a way that t ^ and p^^^i ^ oo as /3 ^ oo. Moreover, 
we will have conditions of the form 

(^1/4^)^0, and r(pi/^£)^oo 



as p ^ oo, where r, v are universal constants. 
Consider now the n-particle Hamiltonian 

n 

^m"a = -iE<,.+^^M, (8) 

where we have introduced the Neumann Laplacian Aq^ ^ of the cube Q^,x = + [—^^j ^^]^ 

and the potential 

Wi,^xixi,...,Xn) = E w{Xi + {fJ. + X)i,Xj + {p + X)i) 
l<i<j<n 



n „ 

-pV] / w{xj + {f^ + X)£,y+{p + X)e) dy 
+\p^ / / w{x + {p + X)£, y + {p + X)£) dx dy. 



3.2. LEMMA (Decoupling of boxes). 

Let EJ^ ^ be the ground state energy of the Hamiltonian H'^ ^ given in considered as a bosonic 
Hamiltonian. The ground state energy Eq of the Hamiltonian H in (01) is then bounded below as 



uj{t)N 



<n<N f"'-" 21 
1 113 



2 ' 2J 

-2/ 



Proof. If ^{xi, . . . ,xn) G L (R ) is a symmetric function. Then 

L0{t)N 



{^,H^) >J2 I (^.^M^) 



2i ' 
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where 

N 



, X]\f)\'^ dxi . . .dxN 



+ 7 / W^^\{xi,...,XN)\'^{xi,...,XN)\'^dxi...dXN- 



The lemma follows since it is clear that {^,H^ \^) > u\ii<in<N ^ □ 

For given ^ the Hamiltonians H"^^ fall in three groups depending on A. The first kind for which 
Qa,^ n = 0. They describe boxes with no background. The optimal energy for these boxes are 
clearly achieved for n = 0. The second kind for which Qa.^* C 17. These Hamiltonians are all 
unitarily equivalent to "fH"^, where 



+ ^ w{xi,Xj) + ^p'^ w{x,y)dxdy, (9) 

l<i<j<n 

where —Ai is the Neumann Laplacian for the cube [— £/2,^/2]^. Finally, there are operators of the 
third kind for which Q^^\ intersects both 0, and its complement. In this case the particles only 
see part of the background. If we artificially add the missing background only the last term in the 
potential W^,a increases. (The first term does not change and the second can only decrease.) In 
fact it will increase by no more than 

^ -"^ ^' ^' w{x,y)dxdy <\p^ [[ \x - y\~^ dx dy < C p^l^ . 



iP j j w(x,y)axay s 2/ 

xG[-t/2,e/2]3 
yel-e/2,i/2]3 

Thus the operator ^ of the third kind are bounded below by an operator which is unitarily 
equivalent to jH^ - Cp^i^. 

We now note that the number of boxes of the third kind is bounded above by C{L/i)'^. The 
total number of boxes of the second or third kind is bounded above by (L + = (1 + L/i)^. 

We have therefore proved the following result. 

3.3. LEMMA (Reduction to one small box). 

The ground state energy Eq of the Hamiltonian H in ^ is bounded below as 
Eo>{l + L/£)S inf inf Spec - C{L/ifpH^ - 

l<n<A' Zl 

where is the Hamiltonian defined in 

In the rest of the paper we shall study the Hamiltonian (M). 
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4 Long and short distance cutoffs in the potential 



The potential in the Hamiltonian @ is ti; given in Our aim in this section to replace w by a 
function that has long and short distance cutoffs. 
We shall replace the function w by 



where 



WrM^^ y) = Xi{x)Vr,R{x - y)xe{y) 

^-\x\/R _ ^-\x\/r 



Vr,R{x)=YR-,{x)-Y,^.{x) 



(10) 



(11) 



Here < r < i? < u;(t)"^^. Note that for x < r then Vr^R{x) r^^ - R^^ and for \x\ > R then 
Vr,R{x) ^ Ixj-ie-l^'l/^. 

In this section we shall bound the effect of replacing w by Wr,R- We shall not fix the cutoffs r 
and R, but rather choose them differently at different stages in the later arguments. 

We first introduce the cutoff R alone, i.e., we bound the effect of replacing w by WR{x,y) = 
Xeix)VR{x — y)Xiiy), where Vr{x) = |x|~^e~l^l/^ = y^-i(x). Thus, since R < uj{t)~^i, the Fourier 
transforms satisfy 

1 1 



Y^/e{k) - VR{k) = in 



> 0. 



(We use the convention that f{k) = f f{x)e~'^^^ dx.) Hence w{x, y)—wji{x, y) = Xi{x) {Y^^^^i — Vr) {x- 
y)Xe{y) defines a positive semi-definite kernel. Note, moreover, that {Y^^/i — Vr) (0) = R~^ —lo/£ < 
R-^ Thus, 



^ w{xi,Xj) - / w{xj,y)dy + \p^ j j w{x,y)dxdy 

I WR{xi,Xj) - j WR{xj,y)dy + j j WR{x,y)dxdy 

\l<i<j<n j=l 

n 

^6{y - Xi) - p 



^5(x -Xi) - p 



{w - Wr){x,y) 



dx dy 



-\jZ^i{x,f {Y^i, - Vr) (0) > -\n {Y^i, - Vn) (0) = -\nR-\ 



(12) 



We now bound the effect of replacing vjr by Wr^R- I.e., we are replacing Yr{x) = |x| 
by (e~l^'l/^ — e~l^'l/''). This will lower the repulsive terms and for the attractive term we get 



n .. n „ 

WR{xj,y)dy > -p"^ Wr,Rixj,y) dy 



> 



r Q-\^-y\/r 

-npsup / Xe{x)- -Xe{y)dy 

X J \x-y\ 

n „ 

-P"^ Wr,R{xj,y)dy -Cnpr'^. 



(13) 
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If we combine the bounds (|T2|) and (|T^) we have the following result. 

4.1. LEMMA (Long and short distance potential cutoffs). 

Consider the Hamiltonian 

j=l V / l<i<j<n 



+ ^p^ J J Wr,R{x,y)dxdy, (14) 

where Wr,R is given in ^1^) and ([TJj with < r < R < uj(t)~^£ and — as before is the Neumann 
Laplacian for the cube [—i/2,i/2\^. Then the Hamiltonian H"^ defined in ^ obeys the lower bound 

> Hl^^ji - \nR-^ - Cmpr^. 

A similar argument gives the following result. 

4.2. LEMMA. With the same notation as above we have for 0<r'<r<R<R'< uj{t)^^i that 

Hlr'^R' > Hl^^R - \nR-^ - Cinpr\ 

Proof. Simply note that Vr'fii {x) — Vr^R{x) = Yj^i-i (x) — Yr-i (x) + Y^-i (x) — Y^i-i (x) and now use 
the same arguments as before. □ 

4.3. COROLLARY (The particle number n cannot be too small). 

There exists a constant C > such that if uj(t)^^ p^^^l > C then HJ^ > if n < Cp£^. 

Proof. Choose R = p^^^^ and r = ^R. Then we may assume that R < uj{t)-'^£ smce uj 



is large. From Lemma 4.1 we see immediately that 

if" > Wr,R{xj,y)dy + ^p^ / / Wr,R{x,y) dx dy - CnpR^ 

.-;::,./.,.:.,.,-..,.//.,..,.«-«. 

The corollary follows since sup^. f Wr^R{x,y) dy < AttR'^ and with the given choice of R and r it is 
ea.y to see that 1 ,)<i.<i,>ci!V^ □ 

5 Bound on the unimportant part of the Hamiltonian 

In this section we shall bound the Hamiltonian H'^^. ^ given in (ll^). We emphasize that we do not 
necessarily have neutrality in the cube, i.e., n and pt^ may be different. We are simply looking for 
a lower bound to i?"^ that holds for all n. The goal is to find a lower bound that will allow us 
to conclude that the optimal n, i.e., the value for which the energy of the Hamiltonian is smallest, 
is indeed close to the neutral value. 
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We shall express the Hamiltonian in second quantized language. This is purely for convenience. 
We stress that we are not in any way changing the model by doing this and the treatment is entirely 
rigorous and could have been done without the use of second quantization. 



Let Up, ip/-K € (N U {0})'^ be an orthonormal basis of eigenfunctions of the Neumann Laplacian 



-A^ such that — A^Up = |p| Up. I.e., 



3 / 

,{xi,X2,X3) = Cpi ^/^Jjcosf 



PjTT{Xj + £/2) 



where the normalization satisfies cq = 1 and in general 1 < Cp < \/8- The function no = ^-3/2 is 
the constant eigenfunction with eigenvalue 0. We note that for p 7^ we have 

{up,-AeUp)>7TH-^. (15) 

We now express the Hamiltonian H^^^ in terms of the creation and annihilation operators 
Qp = a{up) and a* = a{up)* . 
Define 

^pq,tiu = j j Wr,Rix,y)up{x)uq{y)u^{x)u„{y) dxdy. 
We may then express the two-body repulsive potential as 



Wr,R{Xi,Xj) = \Y^ 
l<i<j<n pg,fJ-'^ 



^pq,^v'^p^q^v^^i^ 



where the right hand side is considered restricted to the n-particle subspace. Likewise the Back- 
ground potential can be written 

n 

-pYwr,R{xj,y) dy = -pi^ ^ ^5op,Og«X 

j=l pq 

and the background-background energy 

i/ = 4.' Aw. 

We may therefore write the Hamiltonian as 

Hi,r,R = Yl \P\'^^P^P + kYl ^pg,f^uO'pala^a^ - p£^ ^ ^'op.oqOpOg + kp^^^^oofio- (16) 

p pq,fiu pq 

We also introduce the operators no = agflo n_|_ = Y2p=io- These operators represent the 
number of particles in the condensate state created by Og and the number of particle not in the 
condensate. Note that on the subspace where the total particle number is n, both of these operators 
are non- negative and n_|_ = n — uq. 



Using the bounds on the long and short distance cutoffs in Lemma 4.1 we may immediately 
prove a simple bound on the expectation value of n+. 
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5.1. LEMMA (Simple bound on the number of excited particles). 

There is a constant C > such that if uj{t)~^ p^^^i > C then for any state such that the expectation 
(H^) < 0, the expectation of the number of excited particles satisfies < Cnp^^^^ 



Proof. We simply choose r = R = p ^^'^ in Lemma 4.1, This is allowed since R < uj{t) ^£ is 



ensured from the assumption that uj{t) ^p^^^i is large. We then obtain 

n n 

HI > -i7-iA^j - - Cnpr"^ > ^ -^-f-^A^j - Cnp^/\ 



The bound on follows since the bound on the gap (15) implies that (X]j=i ^ (n+)7r^^ ^. 

□ 

Motivated by Foldy's use of the Bogolubov approximation it is our goal to reduce the Hamiltonian 
^7rR ^° that it has only what we call quadratic terms, i.e., terms which contain precisely two Op 
with p 7^ 0. More precisely, we want to be able to ignore all terms containing the coefficients 

• ""^00,00- 

• ^pQ,qO ~ "^Op.Oq' where p,q ^ 0. These terms are in fact quadratic, but do not appear in the 
Foldy Hamiltonian. We shall prove that they can also be ignored. 

• «5pO,oo = ^opflo = Woo,po = woofip^ where p^O. 

• ^pg,f,o = %o,pg = ^gp,OM = ^0/.,gp> ^here p,q,p^ 0. 

• ''^pq,iiu> where p,q,p,v ^ 0. The sum of all these terms form a non-negative contribution to 
the Hamiltonian and can, when proving a lower bound, either be ignored or used to control 
error terms. 

We shall consider these cases one at a time. 

5.2. LEMMA (Control of terms with i^oo oq). 

The sum of the terms in H^^ ^ containing Wqq qq is equal to 



2'"'oo,oo 



{no - pi^)"^ - no = QQ [n - pfy + {n+y - 2 [n - pl^) n+ - no 



Proof. The terms containing tjjqq qq are 

using the commutation relation [0^,0*] = 6p^q. □ 
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5.3. LEMMA (Control of terms with u5po,qo)- 

The sum of the terms in H"^^ ^ containing Wpo^go or w'op.Og ^^^^ P^Q ¥^ ^ hounded below by 
where = max{t, 0}. 

Proof. The terms containing w^q ^q or Wq^ q^ are 
Note that no commutes with ^ ''^pO,q0^p^q- 



We have that 



w. 



- ,9-3 



'pO,, 



J J Wr,R{x,y) dyup{x)uq{x)dx. 



Hence 



X] "^pO,qO«p«g = £^ Wr,R{x,y)dy {^Up{x)al\ \'^Up{x)a*p^ dx. 



< e ^sup Wr,Rix',y)dy / ^Up{x)al ^Up{x)al dx. 
= sup / Wr,R{x', y) dy ^ 0*0^ = 1'^ sup / Wr,R{x' , y) dyn+. 



Since 



sup J Wr,R{x,y)dy < J Vr,R{y) dy < iirR^ 
we obtain the operator inequahty 

<^WpO qQa*paq < AT:r^B?n+, 

q^tO 

and the lemma follows. □ 

Before treating the last two types of terms we shall need the following result on the structure 
of the coefficients w„„ 

5.4. LEMMA. For all p', q' G (tt/I) (N U {0})^ and a G N there exists J^,^, G M with J^^, = J^y 
such that for all p, g, G {n/i) (N U {0})^ we have 



w 



a 

Moreover we have the operator inequalities 

< X ^PP'.ooOpV = XI ^pO,Op'apV < ^T^i~^R^n+ (18) 
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and 

Proof. The operator A with integral kernel Wr,R{x, y) is a non-negative Hilbert-Schmidt operator on 
L2(m3) with norm less than sup^ Vr,R(k) < A.itE? . Denote the eigenvalues of A by Aq,, q = 1, 2, . . . 
and corresponding orthonormal eigenfunctions by ipa- We may assume that these functions are 
real. The eigenvalues satisfy < Aq < AttB? . We then have 

Wpq^^u = YK / Up{x)u^{x)ipa{x) dx j Uq{y)uy{y)^a{y) dy . 

The identity (|l7|) thus follows with J^^ = A^^ / Up{x)u^{x)ipa{x) dx. 

If P denotes the projection onto the constant functions we may also consider the operator 
(/ — P)A{I — P). Denote its eigenvalues and eigenfunctions by A'^, and cp'^. Then again < A'^, < 
AttR^. Hence we may write 

^pO,Op' =^~'^Y^'a / Up{x)ip'^{x) dx / Up,{y)Lp'^{y)dy. 
Thus, since all Lp'^ are orthogonal to constants we have 

= ^~^YKw1\ Up{x)ip'^{x) dx a*p ] {J2 M^)'PMdxa* 

= r^^^AXK)"^ 

a 

The inequalities ( |l^ ) follow immediately from this. 

The fact that J2pp' m^o ^pm mp'^p%' — ^ follows from the representation (p!7|). Moreover, since 
the kernel 'WR^rix,y) is a continuous function we have that Wr,R{x,x) = X^o '^o'^a(^)^ almost 
all X and hence 

m^O 

We therefore have 



X] ^pm,mp'(^l%' < X] j Up{x)upi{x)Wr,R{x,x)dxalap, 

p,p',mj^O p,p'j^O 

= / Wr,Ji(a;,x) ^'Up(x)a* ^Up(a;)a* dx 

< sup Wr^R{x' ,x') / ^^Up(2;)a* j 'Y^Up{x)a* J = sup 

\P7^0 / \p^0 J 

and the lemma follows since sup^,/ Wr^R^x' , x') < r~^. □ 
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5.5. LEMMA (Control of terms with Wpo,oo)- 

The sum of the terms in H^^. ^ containing w^q^qq, '"5op,oO''"^oo,pO' ^oo,Op' '"^^^^ p ^ Q is, for all 
e > 0, hounded below by 

- e-^ATTr^R^nQn+ - etOoo.ooC^o + 1 - P^^f, (19) 

and by 

- e~'^ATTr^R^non+ - e{^oo,oo(^+ - 1)^- (20) 
Proof The terms containing Wpo^QQ, ^'op,oO'^oo,pO' ^oo,Op are 

= X ^pO,oo ((^0 - P^^)a*paQ + alap{nQ - pi^)) 

= X X "^P^o^o^o (a;ao(no + 1 - pf ) + (no + 1 - /o^^)^^^) • 

a p^O 

In the last term we have used the representation (|l7| ) and the commutation relation [no,oo] = oo- 
For all e > we get that the above expression is bounded below by 



e-^ X E Vh'^oalap, - e ^ {Jmf (^o + 1 - /^^^ 



'3\2 



a p,p't^O 



= -e ^ X ^pcop'^oa^ap, -e^^}oo,oo("^o + l-p^^)^• 
p,pYo 

The bound (^) follows from (|T^). 

The second bound ( pO[ ) follows in the same way if we notice that the terms containing WpQQQ, 
"f^Op.oO'^ocpO' or Wqo.Op ™ay be written as 

X "^pO,oo (('^ - P^^WpaQ + aQQpin - pi^)) 

p^O 

+ X X Voo (opOo(l - n+) + (1 - ^+)aSap) • 

a p^O 

□ 

5.6. LEMMA (Control of terms with Wpg^^o). 

The sum of the terms in HJ^r,R containing Wpq ,^Q, Wpqflrn!''^pO,qm! ^'^ '^Op,qm! p,q,m is 
bounded below by 

-e~^47ri~^R^non+ - en+r'^ -e ^ '^mp',pm'(^*m(J'*p'a^'%, 

p,m,p' ,m'ytO 

for all e > 0. 



EHLJPS-Sept-30-2000—Foldy 



20 



Proof. The terms containing Wp^^^Q, Wp^fim^Wp^ g^, or WQp g^ are 



pqm^O 



+e I ^ ^ Jpj^a^ttp j j ^ ^ jpyy^ttptt^ 



Using that = we may write this as 



qq'y^O p,m,p' ,Tn'j^O 



"^■mp' ,pm' ^m'^p' ^m' '^p 



qq'^0 p,m,p' ,m' ^0 

* 

' mp ,pm' ^ra^m' • 



p,m,m'j^O 



The lemma now foUows from Lemma 5.4. □ 



6 Analyzing the quadratic Hamiltonian 

In this section we consider the main part of the Hamiltonian. This is the "quadratic" Hamiltonian 
considered by Foldy. It consists of the kinetic energy and all the terms with the coefficients Wpg^oQ, 
^oo,pq %o,Og> and WQp gQwith. p,q^O, i.e., 



p 

+ i E ^Pifio (apOoao^g + aSopOgao + apOgOoao + aSaJapag) (21) 
= E \P\^^P^P + E ^P'J.oo («p«g"oao + ^a*pa*gaQaQ + ^a*Qa*Qapag) . 

P pq^O 



In order to compute all the bounds we found it necessary to include the first term in ( |20| ) into the 
"quadratic" Hamiltonian. We therefore define 

= XI \P\'^^P^P + E ^pO,oo ((" - P^^)a*paQ + alap{n - pi^)) 

P p^O 

+ E ^Pifio (opOg^oao + ^Opflgaoao + ^«oaoapag) • (22) 

pqy^O 
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Note that -ffFoidy = Hq in the neutral case n = p£^. Our goal is to give a lower bound on the 
ground state energy of the Hamiltonian Hq. 

For the sake of convenience we first enlarge the one-particle Hilbert space ([— ^/2, £/2]^) . 
In fact, instead of considering the symmetric Fock space over ([— ^/2, ^/2]^) we now consider 
the symmetric Fock space over the one-particle Hilbert space ([— ^/2, £/2]^) © C. Note that the 
larger Fock space of course contains the original Fock space as a subspace. On the larger space we 
have a new pair of creation and annihilation operators that we denote Oq and clq. These operators 
merely create vectors in the C component of ([— ^/2, £/2]^) © C, and so commute with all other 
operators. 

We shall now write 

a = and 2 = | ' ^ ' (23) 

^ \ 5o, ifp = " [ a*o, ifp = 

We now define the Hamiltonian 

= Yl I^'I^^pS + ^pO,oo ((" - P^^KpttQ + aoSp(n - pi^)) 

p p 

+ Yl ^Pifio («P«9«oao + |a;a*aoao + ^a*oa*Qapag) , (24) 

pq 

where we no longer restrict p, q to be different from 0. Note that for all states on the larger Fock 
space for which (Sq^) = we have (Hq) = (Hq). 

For any function ^ L? ([— £/2, £/2]^) we introduce the creation operator 

p 

Note that the sum includes p = 0. the difference from a*{'p) is given by a*[ip) — a*{p) = 
{uQ,p) {al-al). 

Then \fl{ip) ,a* {tp)] = {(p^ip). We have introduced the "dummy" operator Oq in order for this 
relation to hold. One could just as well have stayed in the old space, but then the relation above 
would hold only for functions orthogonal to constants. 

For any A; G M'^ denote Xe,kix) = e'^^^Xe{x) and define the operators 

bl = a* {Xe,k)ao and bf, = a{Xi,k)al 
They satisfy the commutation relations 

[^fe' K'] = OoOo Xt,k') - a{Xe,k)a*{Xe,k') = aQa^xjik' - k) - a{Xe,k)a* {Xe,k') (25) 

We first consider the kinetic energy part of the Hamiltonian. We shall bound it using the double 
commutator bound in Appendix ^ First we need a well known comparisson between the Neumann 
Laplacian and the Laplacian in the whole space. 
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6.1. LEMMA (Neumann resolvent is bigger than free resolvent). 

Let Pg denote the projection in L^(R'^) that projects onto L'^{[—£/2,£/2]^) (identified as a subspace). 
Then if —A denotes the Laplacian on all ofM? and — is the Neumann Laplacian on [— £/2,£/2]^ 
we have the operator inequality 

(_A^ + a)-i > P^(-A + a)-ip^, 

for all a > 0. 

Proof. It is clear that for all / G L'^{m?) 

||P,(-A£ + a)V2p,(_A + a)-V2/f < 

and hence 



||(_A + a)-V2p,(_A, + a)i/2p^/||2 < 
Now simply use this with / = (— + a)^^/^M. 



□ 



6.2. LEMMA (The kinetic energy bound). 

There exists a constant C" > such that ifC't < 1, where t is the parameter used in the definition 
of Xi in Section we have 



|fc|2 + (^t3)-2 



{blbk) dk 



for all states with (aQag) = and particle number equal to n, i.e., \^pa*a^ 

2 



n 



Proof. Let s, with < s < t, be a parameter to be chosen below. Recall that t is the parameter 
used in the definition of Xe in Section ^ Then since x| + ?7| = 1 we have 



Xi 
+ 



Ae + {{£s)-^ 2 



-Ai + {is) 



-2 ' 2 



-A, + {ls)~^ 



A, + {ls)-^ 
-A, + (£s)-2 



(-A£)2 



-A, + (f5)-2 



(-A^)' 
^ ^^ -A, + (£.)^2 ^^ + ^r 

- c{it)-' 



+ 
(-A. 



(-A, 



-A, + 



-A, + (£s)-2 



-A, + 

where the last inequality follows from Lemma B.l in Appendix ^ We can now repeat this calcu- 
lation to get 

(-A,)2 -A. 



-Ai, > Xt 



-A, + (^s)-2 



C{£t) 



-A, + (^s)-2 
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c{et) 



-2 



-A, 



-2 



+ 



-A, 1 



-A, + 



If we therefore use (53) in Lemma p.l| and recall that s < t we arrive at 

(-A,)2 , _A, 



-A£ > X£ 



+ 

Note that for a > we have 

(-A,)2 



A, + (^s)-2 
f-A,)2 



a- 



-A, + {is)-^ 



C{H)- 



C{£ty 



C{lt) 



Xe 



-2 



-A, 



-A, + 



-2„2.-4 



-A, 



-A, + (£s)-2 ^ ^ _A, + (ls)-2 
Thus if we also assume that a < 1 we have 

^-A,)2 



l„2+-4/)-2 



-Ai > (1 - a)Xe- 



-A, + {is)-^' 

Thus if u is a normalized function on ^^(R^) which is orthogonal to constants we have according 
to the bound on the gap (|l5|) that for all < 5 < 1 



f— A )2 

{u,-Atu) > {l-6)il-a){u,Xe^^—^i-^Xeu 



-A, + {lsy^' 

- c{i - 6)a-hh-H-^ + 6TT^r^. 

We choose a = 6 = C st~^'^ for an appropriately large constant C" > and assume that s and t are 
such that 6 is less than 1. Then 

-A,)2 



(n,-A£n) > {l-C'st-y [u,Xi 



;AlX£U 



-Ae + {isY 



;Xeu 



If we now use Lemma p.l| we may write this as 

1 



{u,-Aeu) > {l-C'st~^f ( u,XeA. 



:i-c'st-^? 



i-Af 
^'^^ -A + (^.)-2 ^^^'' 



-A+ (^S)"2 

where in the last inequality we have used that Ax = AiX and xA = xA^. 

We now choose s = and we may then write this inequality in second quantized form as 



M3 |A:|2 + (^t3)- 



{a* {Xe,k)a{Xe,k)) dk 



using that (ogaQ) = 0. Since we consider only states with particle number n the inequality still 
holds if we insert n'^aQUQ as in the statement of the lemma. □ 

With the same notation as in the above lemma we may write 

Wr,Ri^,y) = i2iTy^ / Vr,Rik)Xi,k{x)Xe,k{y)d'k- 
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The last two sums in the Hamiltonian can therefore be written as 



pq 



Note that it is important here that the potential Wr r contains the localization function Xt- 



Thus, since V^,_r(/c) = K-,_r(— A;) and Xt{k) = Xt{—k) we have for states with {aQa^) = that 

(Hq) > I {hQ{k)) dk-Y, ^p,,oo , (26) 

where 



pq 



hqik) 



(1 - C't 



U\2 



{blb^ + bUh-u) 



"2(27r^)3 



6.3. THEOREM (Simple case of Bogolubov's method). 

For arbitrary constants A > B > and k G C uie have the inequality 



Aiblb, + bUb-k) + mb-k + bub-k) + <bl + 6_fe) + K{b^ + b*_,) 
>-l{A- V:4^^)([6„ bl] + b*_,]) ^''^'^ 



A + B 



Proof. We may complete the square 



if 



Aiblb, + b*_,b_,) + Biblbl, + b,b_,) + Kibl + b_,) + + b*_,) 

= D{bl + a6_fc + a){bj. + a^l^ + a) + D{b*_f. + + a){b_i^ + a^l;. + a) 

D{l + a^)=A, 2Da = B, aD{l + a) = K. 



We choose the solution a = A/B — \J A? jB"^ — 1. Hence 



L>(l + Q2 + 2a) ^ + 



(27) 



□ 



Usually when applying Bogolubov's method the commutator is a positive constant. In this 

case the lower bound in the theorem is actually the bottom of the spectrum of the operator. If 
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moreover, A> B the bottom is actually an eigenvalue. In our case the commutator [6^,6^] is not 
a constant, but according to (|25|) we have 



[bk, bl] < j Xiixf dxalao < l^'alao- (28) 
From this and the above theorem we easily conclude the following bound. 

6.4. LEMMA (Lovi^er bound on quadratic Hamiltonian) . 

On the subspace with n particles we have 

where I = ^(27r)-3 f^^ f{k) - {f{kf - gik)^/^ dk with 

1 1 

and 

Proof. We consider a state with (aQOg) = 0. Then (Hq) = (Hq). We shall use (|2^). Note first that 

\ PI I \p,q¥^0 I 



by (ll) and the fact that R < £. We may of course rewrite £ = n^^^£ ^/^( 



nt 



-1/4 



By Theorem 3.3, (27) and (^) we have 



where 



2(27r£)3' •^'^ " ' 4(27r)37n |fc|2 + (^^3)-2• 
Since > Bk we have that 

V;.,i?(A:)(n-p£3)2 



hQik) > -(A - - - ^••^^;ff)3/ • 



Note that 

Vr,Rmn-p£^) 



2(27r)3£6 

1 2 



The lemma now follows from (^) by a simple change of variables in the k integral. □ 
As a consequence we get the following bound for the Foldy Hamiltonian. 
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6.5. COROLLARY (Lower bound on the Foldy Hamiltonian) . 

The Foldy Hamiltonian in MI) satisfies 



HFoldy > -In'^H-^/* - A7:n'/H-'/\ni)-'/\ (29) 

There is cosntant C > such that if p^/'^R > C, p^/'^lt^ > C, and t < then the Foldy 

Hamiltonian satisfies the hound 

HfoMj > i J^lpl'aX-Cn^/^r^/^ (30) 
p 

Proof. Lemma 3.4 holds for all p hence also if we had replaced phy njl& in this case we get (29). 



The integral / satisfies the bound 

/ < ^(27r)-3 / max {g{k), \g{kf{f{k) - g(}i))-^] dk. 

By Corollary O we may assume that n > cp£^. Hence / is bounded by a constant as long as p 



and p^/^£t^ are sufficiently large and t is sufficiently small (which also ensures that 7 is close to 1). 
Note that we do not have to make any assumptions on r. Moreover, if this is true we also have that 
n£ > cp£'^ is large and hence {n£)~^ is small. This would give the bound in the corollary except 
for the first positive term. The above argument, however, also holds (with different constants) if 
we replace the kinetic energy in the Foldy Hamiltonian by ^ (7"^ — ^) Ipp^p'^p (assuming that 
7 < 2). This proves the corollary. □ 

Note that if 

^1/4^-3/4^ ^ cx), n^/^r^/V ^ 0, ^1/4^1/4^3 ^ and t ^ (31) 
it follows by dominated convergence that / converges to 

i(27r)-3 / iTT\k\~^ + i|A:|2 - {i47T\k\~'' + l\k\Y - (4^1^1-2)^)'/' dk 

= (2/7r)3/4 1 + _ ^2 (^4 ^ dx = - (^^^A, 

where A was given in (^. Thus if we can show that n ~ p£^ we see that the term —In^/'^£~^/^ ~ 
—Ip^^^n agrees with Foldy's calculation for the little box of size £. 

Our task is now to show that indeed n ~ p£^, i.e., that we have approximate neutrality in each 
little box and that the term above containing the integral / is indeed the leading term. 

7 Simple bounds on n and n+. 



The Lemmas 1.1,5^^,5.3,5.5, and 5.6 together with Lemma or Corollary |6.5| control all terms in 



the Hamiltonian H? except the positive term 



2 ^ ^ "^mp' ,pm''^m^p'^m''^p• 
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If we use (^) in Corollary together with the other bounds we obtain the following bound if 
p^^^R and p^/^lt^ are sufficiently large and t is sufficiently small 



"2^00,00 



(no - ) - no 
-47r[p - n£-3]^n+ii2 - Annlr^R^ 
-e'^STTl^^R^non^ - e'Woo,oo("'0 + 1 - P^^)^ 
-en+r~^ + (l - f^) X] ^mp'.pm'amap'Om'ap- 

The assumptions on p^^'^R, p^/^it^, and t are needed in order to bound the integral / above by 
a constant. If we choose e = 1/4, use i^oooo — ^vri?^^"^ and ignore the last positive term in the 
bound above we arrive at 

m > l^\p\'a;a^-Cn'^'i-^/^-^nR-'-Cnpr^ + lw,,^oo{'^o-pi'^^ 
p 

-AttIp - n£-^]+n+R^ - Ai:n\r^R^ 

-32Tr£-^R'^non+ - AttR^T^ (no - ^P^^ + 3) - iKr'^ 

> 3 E \P\^''>P - C™'^'^"'^' - h^R~' - Cnpr^ + i{2oo,oo (^o - P^')' 
p 

-A87ri-^R^nn+ - AttR^I'^ (no + 3) - 3^+r"\ (32) 
where in the last inequality we have used that pl"^ < 2n, uq < n and n+ < n. 
7.1. LEMMA (Simple bound on n). 

Lei Lj(t) 6e t/ie function described in Lemma 3.1. There is a constant C > such that if {p^^^i)t^ > 



C and {p^/^tjp 1/12^ Q^jid uj{t){p^/^£) ^ are smaller than C ^ then for any state with (H^) < 
we have C~^p£'^ < n < Cp£^. 



Proof. The lower bound follows from Corollary [4.3| . To prove the upper bound on n we choose 
R = uj{t)^^£ (the maximally allowed value) and r = bio{t)^^£, where we shall choose b sufficiently 
small, in particular b < 1/2. We then have that p^^^R = uj{t)~^ p^^^£ is large. Moreover Woqoo — 
CR?'£~^ = Cuj{t)^'^£^^ for some constant C > and we get from ( |3^ ) and Lemma |5.1| that 

{H^) > £-^[-Cn^/'^£^/^ - ^nuj{t) - Cb^uj{t)-^n^ + CLo{t)-mno) - p£^)^ 
-A8Truj{t)-^p-^/^{£p^/yn^ - 4Truj{t)-^ + 3) " inb-^Lo{t)], 

where we have again used that cp£^ < n, no < n and n^ < n. Note that 

n^/^/^ < Cu;{tr\\p'/^£)-^p-'/Utf 

and nLo{t) < Cuj{t)^'^n^ p^^uj{t)'^ From Lemma |5l] we know that (no) > n(l — Cp~^/^(£p^/^)^). By 
choosing b small enough we see immediately that n < C p£^. □ 
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Using this result as an input in (32) we can get a better bound on n than above and a better bound 



on (n^) than given in Lemma 5A. In particular, the next lemma in fact implies that we have near 



neutrality, i.e., that n is nearly pi^. 

7.2. LEMMA (Improved bounds on n and 

There exists a constant C > such that if {p^/'^i)t^ > C and (p^/^^)/?"^/^^, t, and uj{t){p^/'^iy^ 
are smaller than then for any state with (H^) < we have (X^„ IppflpOp) 
and ^ 

{n+)<Cnp-^/\p^/Hf and (^^^) < C p-^'\p^lHf . 

For any other state with {Hf^., j^,)' < we have the same bound on {n+Y if r' < p-3/8(-^i/4£)i/2 
and R' > a{p^/^t)^'^i where a > is an appropriate constant. 



Proof. Inserting the bound n < Cpi"^ into ( |3^ ) gives 

p 

-CR^pn+ - CR^ {p + le-^) - ln+r-\ 

We now choose r = p-3/8(^i/4^)i/2 ^ ^ a{p^/^e)-^£, where we shall choose a below, in- 
dependently of p, p^^'^i, and t. Note that since uj{t){p^^^£)~'^ is small we may assume that 
R < uj(t)^^£ as required and since is smah we may assume that r < R. Moreover 

r-i = /,-i/8(pi/4^)3/2^-2 j^2p ^ a^{p^/H)-H'^p = a^r^. Hence, since Yip IPpa^^ > T^'^r'^n+ 
(see 15), we have 

_(^ + op'^/H^ip^'H) - ca^p'>/H\p^'Hr\i + {p^Hr^'p-"^). 

By choosing a appropriately (independently of p, p^^^l., and t) we immediately get the bound on 
( IppCpOp) and the bound ^~^(n+) < C p^/'^i^{p^^'^l), which implies the stated bound on (n+). 
The bound on (n — pi^)'^ {p£^)~'^ follows since we also have Wqq qq 
and 

WooMi'^o-pi'f) > CR'i-mno) - p£^f 

> Ca\p^lHrH^ (n - pe - nCp-^l\£p'ly)\ 

where we have used the bound on which we have just proved. 

The case when {Hf^, j^,)' < follows in the same way because we may everywhere replace 



by H^^, ^, and use Lemma 4.2 instead of Lemma 4.1. Note that in this case we already know the 



bound on n since we still assume the existence of the state such that {H^) < 0. □ 
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8 Localization of n_ 



Note that Lemma 7.2 may be interpreted as saying that we have neutrahty and condensation, in 
the sense that is a smah fraction of n, in each Httle box. Although this bound on (n+) is 

sufficient for our purposes we still need to know that (n^) ~ We shall however not prove 

this for a general state with negative energy. Instead we shall show that we may change the ground 
state, without changing its energy expectation significantly, in such a way that the possible n+ 
values are bounded by Cnp-'^/^{p^/^£f . To do this we shall use the method of localizing large 
matrices in Lemma |A.1| of Appendix ^ 

We begin with any normalized n-particle wavefunction ^ of the operator H^. Since ^ is an 
n-particle wave function we may write ^ = X]m=o '^rn^m, where for all m = 1, 2, . . . , n, is a 
normalized eigenfunctions of n+ with eigenvalue m. We may now consider the (n + 1) x (n + 1) 
Hermitean matrix A with matrix elements .4mm' = \^m,H'^r R'^'m 



We shall use Lemma A.l for this matrix and the vector ^/^ = (cq, . . . , c^). We shall choose M in 



Lemma A.l to be of the order of the upper bound on l(n.+) derived in Lemma 7.2, e.g., M is the 



integer part of np ^/^{p^/'^i)'^ . Recall that with the assumption in Lemma 7.2 we have M ^ 1. 



With the notation in Lemma A.l we have A = {tp,A^) = , H^^. j^"^). Note also that because of 



the structure of H^^ we have, again with the notation in Lemma [A.lJ , that = if /c > 3. We 



conclude from Lemma A.l that there exists a normalized wavefunction ^ with the property that 



the corresponding values belong to an interval of length M and such that 

{^,Hl,^j,^) > (^,Hl,^j,^) - CM-\\d,\ + \d2\). 



We shall discuss di,d2, which depend on in detail below, but first we give the result on the 
localization of n_(_ that we shall use. 

8.1. LEMMA (Localization of 

There is a constant C > with the following property. If {p^/^i)t^ > C and t, and 

a;(t)(pi/^^)-i are less than C'^ and r < p3/8(pi/4^)i/2^ ^ > C{p^/^£)-'^l , and ^! is a normalized 
wavefunction such that 

{^,Hl,^R^) < and < -C{np-^/\p^/Hf)-Wdi\ + jdal) (33) 

then there exists a normalized wave function ^, which is a linear combination of eigenfunctions of 
n+ with eigenvalues less than Cnp^^^^{p^^^i)^ only, such that 

(*,/?," pM/) > - C{np-y\p'/'if r\\d,\ + m. (34) 



Here di and d2, depending on , are given as explained in Lemma A.l. 



Proof. As explained above we choose M to be of order np ^/^(yO^/^£)^. We then choose as 
explained above. Then (p3) holds. We also know that the possible n+ values of ^ range in an 
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interval of length M. We do not know however, where this interval is located. The assumption 
(^) will allow us to say more about the location of the interval. 

In fact, it follows from (p3|), ( p^ that Hf^j^^^ < 0. It is then a consequence of Lemma [7^ 

that (^$,n+*) < Cnp-^/'^ip^/'^if. This of course establishes that the allowed n+ values are less 
than C'np^'^/^{p^/^ef for some constant C" > 0. □ 

Our final task in this section is to bound di and d2- We have that di = , H^^ j^(l)'tlj), where 



w. 



pq,fj.u 



H"^^ Plil) is the part of the Hamiltonian H'^^j^ containing all the terms with the coefficents 
for which precisely one or three indices are 0. These are the terms bounded in Lemmas and 5^ . 
These lemmas are stated as one-sided bounds. It is clear from the proof that they could have been 
stated as two sided bounds. Alternatively we may observe that H'^^piil) is unitarily equivalent to 
—H^^j^{l). This follows by applying the unitary transform which maps all operators a* and 



with p 7^ to —a* and —a^. From Lemmas |5.5| and |5.6| we therefore immediately get the following 
bound on di. 

8.2. LEMMA (Control of di). 

With the notation above we have for all e > 

\di\ < £-^8TTr^R^{^,non+^)+e{^,{n+r~^ + WQo^Qo{no + l- pi^f)^) 

\ p,m,p',m'^0 

Likewise, we have that d2 = (^', -ff"^ ^(2)^/>), where H"^j^(2) is the part of the Hamiltonian 
^irR containing all the terms with precisely two Oq or two Oq. i.e., these are the terms in the 
Foldy Hamiltonian, which do not commute with n+. 

8.3. LEMMA (Control of da)- 

There exists a constant C > such that if {p^/'^i)t^ > C and t, and w(t)(p^/^^)"^ 

are less than C^^ and is a wave function with {'^,H^'$) < then with the notation above we 
have 

\d2\ < Cp^/H^{p^/^£) + Airr^R^ (^,n+no^) • 

Proof. If we replace all the operators a* and with p 7^ in the Foldy Hamiltonian by —ia* and 
iUp we get a unitarily equivalent operator. This operator however differs from the Hamiltonian 
-f^Foidy only by a change of sign on the part that we denoted H"^^ ^{2). Since both operators satisfy 
the bound in Corollary ^?5| we conclude that 



< ^ 



P pqi^O 
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Note that both sums above define positive operators. This is trivial for the first sum. For the 



second it fohows from (18) in Lemma 5A since Oq^q commutes with aU a* and with p ^ 0. The 
lemma now follows from (18) and from Lemma 7.2. □ 



9 Proof of Foldy's law 

We first prove Foldy's law in a small cube. Let ^ be a normalized n-particle wave function. We 
shall prove that with an appropriate choice of £ 



^1/4 



.1/4 



(35) 



where A is given in (|6|). Note that ^ < 0. It then follows from Lemma 3^ that 

Eo>{l + L/ifj iff Api' (pV4 + , (^1/4)) _ ciL/ifpH' - 

Thus, since = pL"^ we have 

^lim § > 7 iff A (pV4 + , (,1/4)) _ cp^^u^it) [py^y . 

Foldy's law (|5|) follows since we shall choose (see below) t and £ in such a way that as p ^ oo we 
have t — > and hence 7 — > 1 and uj{t){p^/'^£)~^ — > (see condition (41) below). 

It remains to prove (^Sl). First we fix the long and short distance potential cutoffs 



R = io{t)-^i, and r = />-3/8(/5V4^)-V2 



(36) 



We may of course assume that {^,HJ^^) < 0. Thus n satisfies the bound in Lemma 7.2. We 
proceed in two steps. In Lemma |9.1| Foldy's law in the small boxes is proved under the restrictive 
assumption given in (|3^) below. Finally, in Theorem |9.2| Foldy's law in the small boxes is proved 
by considering the alternative case that (|37| ) fails. Let us note that, logically speaking, this could 
have been done in the reverse order. I.e., we could, instead, have begun with the case that ( |37| ) 



fails. At the end of the section we combine Theorem 9.2 with Lemma |3.3| to show that Foldy's law 
in the small box implies Foldy's law Theorem |1.1| . 



At the end of this section we show how to choose £ and t so that Theorem 9.2 implies (|35| ) and 
hence Theorem |1.1| , as explained above. 



9.1. LEMMA (Foldy's law for H^: restricted version). Let R andr be given by i\3q). There 
exists a constant C > such that if {p^^^£)t^ > C and {p'^l'^£)p~^l^'^ , t, and uj{t){p^^^£)~^ are less 
than C^^ then, whenever 



n£--'R-' (^,n+^') 



(37) 
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we have that 



with I as in Lemma \6.A 



Proof. We assume {^,H'^^) < 0. We proceed as in the beginning of Sect. but we now use (|29| ) 
of Corollary instead of (|30|) . We then get 



+5"^oo,oo {no-pff-no 
-en+r-^ + - e) ^ 

p,m,p' ,m'^0 



^mp' ,pm' ^p' ^m' ^p " 



If we now use the assumption (|37|) and the facts that n+ < n, no < n, and w^q qo — 47ri?^£ we 
see with appropriate choices of e and C that 

H"^ > -In^'^r^'^-ATnt'/^r^'^{ni)-^/^-\nR-^-Cnpr'^-CR^r^{n + l) 
-CR^£-^\n - pf\{n+ + 1) - Cn+r-\ 



If we finally insert the choices of R and r and use Lemma we arrive at the bound in the 
lemma. □ 

9.2. THEOREM (Foldy's law for H'^). 

There exists a C > such that if {p^/^£)t^ > C and {p^/'^£)p-^/'^'^ , t, and u;{t){p^/'^£)-'^ are less 
than then for any normalized n-particle wave function ^ we have 

{^,Hf^) > -In^'lH-^l^ -Cp^lH^ {u:{t){p^lH)-^ ^u^{t)-^p-'^l^\p^lHf''l^ 



where I is defined in Lemma 6.4 with r and R as in (p^ 



Proof. According to Lemma 9.1 we may assume that 

n£-^R^ (39) 



where C is at least as big as the constant in Lemma 9.1. We still assume that (^', H^^) < 0. 
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We begin by bounding di and ^2 using Lemmas 3.2 and 8.3. We have from Lemmas 7.2 and 
8.5 that 



\d2\ < Cp'>iH\p"H) + crMtr^n''p-^'\p'/Hf 
< C{np-^i\p'iHffp''iH^^{tY''{p'lH)-\ 



In order to bound d\ we shah use (|3^). Together with Lemma B.2 this gives (choosing e = 1/2 say) 



Inserting the choices for r and R and using Lemma |7.2| gives 

< C{np-^l\p^lHffp'^!H'' (a;(t)-2(pi/^£)-^ + p-i/8(pV^^)-i^/2 



where we have also used that we may assume that p-^l^{p^l^iY'^l'^ is smah. The assumption (|3|) 
now reads 

If this is not satisfied we see immediately that the bound (|38|) holds. 



Thus from Lemma S.l it follows that we can find a normalized n-particle wavefunction ^ with 

%,n^^^<Cnp-^l^{p^^^^f and ,n\^^ <Cn^ p-^l'^{p^l^lf (40) 



such that 



ii7,,K^) > , H^,r,R^) - Cp^l'e (^(t)-2(/,V4^)-7 + ^-l/8(^l/4^)-17/2^) 



In order to analyze \^ ^W^,^ we proceed as in the beginning of Sect. |^ This time we use 
Lemmas 1.1 5^.5.3,5^, and |5.6| together with Lemma 3^ instead of Corollary We obtain 



^7,r,R ^ ^^00,00 



(n - pi'f + (n+)^ - 2 (n - p^^) n+ - no 

p,m,p' ,m' ^0 

_ 1 (n - ^^^^^^ _ 4vrn5/4^-3/4(„^)-l/4 _ J„5/4^.3/4_ 

This time we shall however not choose e small, but rather big. Note that since Wr,R{x,y) < r^^ we 
have ^mp' pm''^*mK)'^m'^p — r~^n^{n^ — 1), which follows immediately from 



p,m,p' ,in' ^0 



E 

p,m,p' ,m' j^O 



^mp' ,pm' ^p' ^m' ^p 



{y)(^m%\ Y Um{x)up{y)a.^apdxdy. 



,p,rn^O 



p,m^O 
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We therefore have 



no 



n Hj 



en:r 



^2 —1 

en , r 



If we now insert the choices of r and R, take the expectation in the state given by ^, and use 
and the bound on n from Lemma 17.21 we arrive at 



> 



_J„5/4^-3/4 _ ^^5/4^3 ^(pl/4£)-l + ^(t)-2(pl/4^)-l 



If we now choose e = uj{t) ^ p^^^^(p^/^i) we arrive at (38). 



□ 



Completion of the proof of Foldy's law, Theorem Ll_. We have accumulated various errors and we 
want to show that they can all be made small. There are basically two parameters that can 
be adjusted, ^ and t. Instead of £ it is convenient to use X = p^l'^l. We shah choose X as a 
function of p such that X — > oo as p — > oo. From Lemma 7.1 we know that for some fixed 

I) we have that 



C > C ^ p(^ < n < Cpi^. Hence according to (31) with r and R given in 



as p ^ oo if 



vo{t)-^X 
p^l'X 
t^X 



t 



oo 

CXD 
CXD 

0. 



(41) 
(42) 
(43) 
(44) 

(45) 

hold. From Lemma |7.2| , for which the hypotheses are now automatically satisfied, we have that 
n = pi^{\ + 0{p~^/^X^/'^) and from ( pSj ) we see that n is pl^ to leading order. 

With these conditions we find that the first term on the right side of (^) is, in the limit p 
exactly Foldy's law. The conditions that the other terms in (|38|) are of lower order are 



The hypotheses of Theorem 9^ are valid if (|4l|),(|4^), (^4|), and 



oo, 



{X/u{t)f'^''p-^I^^^X ^ 



(46) 
(47) 



together with (|4l|). 

It remains to show that we can satisfy the conditions (pT|-47). Condition ([42| ) is trivially satisfied 
since both p and X tend to infinity. Since oj{t) ~ for small t we see that (43) is implied by 
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(^). Condition ( [45D is implied by (47), which is in turn implied by ( pl[ ) and (|4^). The remaining 
two conditions (^) and (p^ ) are easily satisfied by an approriate choice of X and t as functions for 
p with X ^ oo and t ^ as p ^ cxd. In fact, we simply need pi/ii6^-i6/29 ^ x ^ t^^. 



The bound ( pq ) has now been established. Hence Foldy's law Theorem 1.1 follows as discussed 
in the beginning of the section. □ 



A Appendix: Localization of large matrices 



The following theorem allows us to reduce a big Hermitean matrix, A, to a smaller principal 
submatrix without changing the lowest eigenvalue very much. ( The k^^ supra- (resp. infra-) 
diagonal of a matrix A is the submatrix consisting of all elements ai^i^k (resp. ai^k,i)- ) 

A.l. THEOREM (Localization of large matrices). Suppose that A is an N x N Hermitean 
matrix and let A^ , with k = 0, 1,...,A^ — 1, denote the matrix consisting of the k^^ supra- and 
infra- diagonal of A. Let t/j G be a normalized vector and set = {ip,A^il)) and X = {t/^jAip) = 
Ylk=o '^k- (V' need not be an eigenvector of A.) 

Choose some positive integer M < N . Then, with M fixed, there is some n £ [0, — M] and 
some normalized vector (j) G with the property that (pj = unless n-\-l<j<n-\-M (i.e., (j) 
has length M ) and such that 

„ M-l N-l 

{ct,,AcP)<X + j^Yl k^\dk\ + CY, 141 , (48) 

k=l k=M 

where C > is a universal constant. (Note that the first sum starts with k = 1.) 

Proof. It is convenient to extend the matrix Aij to all — oo < i,j < +oo by defining Aij = unless 
^ ^ i,j ^ N. Similarly, we extend the vector ^/^ and we define the numbers dk and the matrix A^ 
to be zero when k ^ [0, N — 1]. We shall give the construction for M odd, the M even case being 
similar. 

For s G Z set /(s) = Am[M + 1 - 2|s|] if 2|s| < M and f{s) = otherwise. Thus, /(s) / for 
precisely M values of s. Also, f{s) = /(— s). Am is chosen so that J2s /(^)^ ~ 

For each m £ Z define the vector ^^'^^ by = f{j — m)iljj. We then define K^'^^ = 

((/)(™),^0(™)) -{\ + cj) ((/)("),(/)("")). (The number a will be chosen later.) After this, we define 
K = T.m Using the fact that Y.s fi^f = 1' ^^ve that 

^(0(™),^,/.(-)) = ^j;(0(-),^V"))=5]^/(s)/(A: + s)(V,^''V') 

m m k=0 s k 

s k=0 

and 

A = A 5^ (</>(-), 0(-)) = J]/(.)2(V,^'=V) = EE/(^)'4 (49) 

m s k=0 s k 
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Hence 

iV-l 

i^ = ^K(-) = -a-^47A: (50) 

m k=l 

with 

Let us choose a = -Y.k=i dklk- Then, Em^^'"'' = 0- Recahing that not aU of the (/>(™) 
equal zero, we conclude that there is at least one value of m such that (i) (j)^'^'^ ^ and (ii) 
((/)(™),^(^("^)) < (A + a)((A(™),(/'(™)). 

This concludes the proof of ( ^8| ) except for showing that < C'^s^jgr foi" all M and k. This 
is evident from the easily computable large M asymptotics in ([5l|). 

□ 



B Appendix: A double commutator bound 



B.l. LEMMA. Let — Ajv he the Neumann Laplacian of some bounded open set O. Given 9 S 
C°°(P) with supp \V9\ C O satisfying \\die\\ < Ct-^ , WdidjOW < Ct''^ , \\didjdke\\ < Ct'^ , for some 
< i and all i,j, k = 1,2, 3. Then for all s > we have the operator inequality 



(-A 



N) 



An + s' 



> -ct- 



-A 



N 



Atv + s- 



Cs't 



2+-4 



We also have the norm bound 



-A 



TV 



Atv + s- 



Proof. We calculate the commutator 



(-A 



Atv + s- 



1 



+ 



-Aiv + s 

-Ajv 
Aat + s- 



<c{sH-^ + sH-^) 



-Atv,! 

-A^,6 



-2 



Aat + s 



-2 



-A 



TV J 



Likewise we calculate the double commutator 



-A 



-An + s-^' 

+ [[-AN,e] 9] 



2s" 



-Atv + s 
-Atv , -A;v 



[[-A^,^] 9] 



-A 



N 



-A 



AT + S 



-2 



Aat + S' 



+ 



Aat + s- 



-A 



N, 



1 



Aiv + s- 



Atv + s 



-2 



Atv + s- 



Note that [[— A^v,^] 9] = —2 (V0)^ and thus the first term above is positive. 
We claim that 

[-Aiv, 9] [9, - Ajv] < -Ct^'^^N + Ct'^. 



(52) 



(53) 



(54) 



(55) 
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To see this we simply calculate 

3 

[-A^, e] [e, -Ajv] = -J2 {^di{die){dje)dj + {dfd){d]e) + 2{die){did]B)) 

The last two terms are bounded by Ct~^. For the first term we have by the Cauchy-Schwarz 
inequality for operators, BA* + AB* < e~^AA* + eBB*, for all e > 0, that 

3 3 3 

and this is bounded above by — St^^A^r and we get (p5[). Inserting ( [55| ) into (^), recalling that 
the first term is positive, we obtain 

\2 



-An + s- 



> -2{V0f 
-Ct 



-A 



N 



-A 



N 



2 -^Af ri„2^-4 



An + s- 



CsH-\ 



Again using the Cauchy-Schwarz inequality, we have 



2{V9f 



< 2r 



An + s~ 



+ 2- 



An + s- 



;{V9Y 



A \ 1/2 



-An + s-2 



-A 



N 



-A7V + S-2 



1/2 



+ 2t' 



-A 



N 



-An + s-2 



< ct- 



-A 



N 



An + s- 



and (ll) follows. 

The bound (^) is proved in the same way. Indeed, 



-A 



N 



An + s- 



1 



+ 2s 



-2 



1 

-A/v + s-2 



Ajv + s- 



:[[-^N,e], 



-An,0] 



1 



AT -I- 6' " — Aat + S^ 

2 



A^ + s-2 
J [0 - An] 



-An + s-^-' 



and (H) follows from [[-Ajv, 6*] 9] = -2 {V9y and (| 



□ 
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